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ABSTRACT 



Generalising recent studies on the sawtooth lattice, a two-spin variant of the model is 
considered. Numerical studies of the energy spectra and the relevant spin correlations in 
the problem are presented. Perturbation theory analysis of the model explaining some of 
the features of the numerical data is put forward and the spin wave spectra of the model 
corresponding to different phases are investigated. 
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Chapter 1 
Introduction 



There has been a lot of recent interest ([4], [5] to [8]) in one dimensional and quasi- 
one-dimensional quantum spin systems having two different spins in the unit cell with 
antiferromagnetic couplings. Depending on the prescence or the abscence of frustration 
and its strength when it is present these systems exhibit a rich variety of phases in 
the ground state. Quantum ferrimagnets for example, are one such class of systems 
where the system has a finite magnetic moment in the ground state. Chemists have been 
successful in synthesising families of organo- metallic compounds ([9] to [12]) which provide 
experimental realisations for some such systems. 

In this report we present the study of a mixed spin variant of the Sawtooth Lattice. 
Recent studies ([1], [2], [3]) on this model have concentrated on systems where all the 
spins on the lattice are the same. The compound Delafossite {YCu02.5) provides an 
experimental realisation of such a model with the copper ions forming a lattice of spin-| 
sites. In an attempt to generalise such studies we considered a two spin variant of the 
above model and studied it numerically and analytically. In this chapter the model and its 
Hamiltonian are introduced and some of the classically expected properties of the ground 
state are discussed. In subsequent chapters the numerical data and the analytical results 
obtained have been presented. 
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1.1 The model and its classical ground states 



The model under consideration is described by the Hamiltonian, 



n n 

Or equivalently, 

H — ~^ln • ~^l{n+l) + X^("^2n " "^In + 2n ' "^l(n+l)) 



(1.1) 



(1.2) 



where S = J1/J2 and J2 has been set to 1. Thus all energies in the problem are measured 
in units of J2. Here Su denotes a spin-1 site and 5*21 denotes a spin-|. 
Schematically the model looks hke, 



spin 1/2' s-> 



spin I's -> 




Classically, the ground state of this system is characterised by a planar or a collincar 
arrangement of spin vectors. Which of these arrangements is the ground state depends 
on the relative strengths of the two couplings Ji and J2. 

When J25'2 > 2Ji5'i the classical ground state is characterised by a coUinear ar- 
rangement of spins and this phase is called the ferrimagnetic state. Schematically the 
ferrimagnetic phase looks as follows: 



When 725*2 < 2JiS'i, the classical ground state is characterised by a planar arrangement 
of spins in each triangle and this phase is called the spiral/canted phase. Schematically 
this phase looks as follows: 




where cos^ = J2S2/2J1S1. Because of the freedom of choosing the direction of spins 
even when the spin vectors in each triangle are constrained to be on a plane, the classical 
spiral/canted phase has an infinite amount of degeneracy. We note that because of the 
above reason though the spin vectors in a particular triangle have to be in a plane, all 
the spin vectors need not lie in the same plane. 

For our simulations and the analytical results that follow, the parameter in the 
problem is S. Numerically we have studied the energy spectra and the spin correlations as 
a function of S. These results are presented in the next chapter. And perturbation theory 
calculations to compute the effective Hamiltonian between the spin-|'s in the large S limit 
and the spin wave spectra obtained for the above two phases are presented in chapters 3 
and 4 respectively. 
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Chapter 2 

Exact Diagonalization and related 
results 

Numerical studies of the model have been done using the Exact Diagonalisation method. 
By calculating the required eigenvalues and eigenvectors using the Lanczos algorithm the 
following quantities were calculated: 

• The variation of the ground state and the first excited state energies with 6. 

• The correlations between the spins in the ground state. 

• The effective Hamiltonian governing the spin-|'s when the coupling between the 
spin-l's is much stronger than J2, the spin-1 - spin-i coupling . 

In this chapter we begin with a brief introduction to the Lanczos algorithm and 
the version of it which has been used. Following that the numerical results of the first 
two categories above are presented. The effective Hamiltonian calculations, being semi- 
analytical in nature are presented in a subsequent chapter. 
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2.1 The Lanczos Algorithm 



The Lanczos Algorithm is a widely used method for finding a few eigenvalues of a large 
symmetric matrix. Since the matrices that one deals with in quantum spin systems are 
usually symmetric (or Hermitian,whose eigenvalue problem can be formulated as one of 
a symmetric matrix double the size), this method is commonly used to find the lowest 
few eigenvalues and eigenvectors corresponding to the ground state and the lowest excited 
states. 

The basic content of the version of the Lanczos procedure used for our simulations 
is the following recursion relation. 

Pi+iVi+i = Avi - aiVi - PiVi-i (2.1) 

for i=l,2 

where A is the matrix whose eigenvalues we want to calculate and ctj = v^Avi and = 
vf+iAvi. Pi is taken to be and vi is chosen to be a random vector normalised to 
unity. For any i — m the a symmetric tridiagonal matrix is defined whose diagonal 
elements are ctj and the off-diagonal elements are (5j (j=2,m). It can be proved (for 

infinite precision arithmetic) that if Ai >= A2 >= A3 >= A^ are eigenvalues of 

and Ai >= A2 >= A3 >= A^ are the m largest eigenvalues of A, then the sequence 

of Aj's converges to the sequence of Aj's as m is incremented. 

Again, let Vm be the nxm (where n is the order of A) whose ith column is Vi. Then 
if Xi is the eigenvector of Aj and Xj, is the eigenvector of Aj then V^Xj — > Xi as m is 
incremented. The vectors VmXi are called Ritz vectors. 

The Lanczos vectors Vi which are generated by the above recursion are an orthonor- 
mal set. The proof of the two claims made above hinges on this orthonormality of the 
Lanczos vectors. In reality when these vectors are generated numerically, finite preci- 
sion effects enter and the set of vectors generated are not strictly orthogonal. The loss 
of orthogonality of the Lanczos vectors affects the computations primarily in two ways. 
Eigenvalues which are simple appear as multiple eigenvalues of the system and more im- 
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portantly spurious eigenvalues which are not eigenvalues of A at all appear as eigenvalues 



Various methods have been devised to overcome the difficulties posed by this loss 
of orthogonahty. One way is to resort to reorthogonahsation of the Lanczos vectors. 

We don't use this method in our computations. We instead use the identification test 
developed by Cullum and Willoughby ([14]) to explicitly identify the eigenvalues which 
are spurious and discard them as they are detected. 

As already mentioned in finite precision calculation the appearance of an eigenvalue 
as a multiple eigenvalue of the tridiagonal matrix is no guarantee of that eigenvalue 
being a true multiple eigenvalue of A. But this difficulty can be overcome by looking at 
the corresponding Ritz vectors. If Ritz vectors are calculated for a large enough Tm for 
which a particular eigenvalue is duplicate then for a true multiple eigenvalue two linearly 
independent Ritz vectors can be generated using appropriate T^'s of different sizes. But 
if the eigenvahic is simple any two Ritz vectors of the same eigenvalue will essentially 
be the same (upto a sign). In this way we can determine the degeneracy of eigenvalues 
by computing more and more Ritz vectors and checking for linear independence. In our 
program we employ this method to determine the degeneracy of an eigenvalue. 



2.2 Results of Exact Diagonalisation 



All the computations have been done by generating the Hamiltonian in the total 
basis. The variation of various energies and correlations with the ratio of interaction 
strength S — J1/J2 has been studied. In all the results reported 6 varies from 0.1 to 
2. The correlations have been calculated for a system size of 10 triangles and all the 
other graphs are for a system size of 8 triangles. All through the computations periodic 
boundary conditions have been used (the spin-l's are joined in a ring). For all the reported 
data the accuracy measured by \{AX - AXy {AX - AX)\ for an eigenvalue A and its 
corresponding eigenvector X (normalised to unity) is ~ 10^^^ or less. 



11 



2.2.1 Variation of ground state energy and first excited state 
energy with S 



Variation of ground state energy with 8 

-8 1 \ \ \ \ \ \ \ — 
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Variation of first excited state energy witti 8 
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2.2.2 Variation of total spin of the ground and first excited 
states with 6 
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2.2.3 Variation of energy gap to the 1st excited state with 6 

Variation of energy gap to 1st excited state with 8 
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Significant aspects of the above results are the following: 



• For a small 5 this sawtooth lattice can be approximated by an alternating spin- 
l/spin-| chain. This system has a ferrimagnetic ground state of total spin N{Si — S2) 
and the first excited state is a state of total spin one less than that of the ground 
state ([4]). One can see clearly from the figures that for small S this is indeed the 
case for this system. 

• There is a sudden change in the total spin of both the ground state and the first 
excited state at around 5 — 0.25. This is the point where we expect the transition 
from ferrimagnetic to the spiral phase from the classical analysis. We have analysed 
this particular region closely by studying the total spin behaviour of the ground 
state at a number of closely spaced points from 5 — 0.2 to 5 — 0.35. Numerically 
we have found that for the ground state the spin drops to zero ai S — 0.265. 
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• There are two values of 6 where the system seems to be gapless. The first is near 
S = 0.5 and the other is near 6 = 1.0. This actually divides the phase diagram into 
three phases as opposed to the two phases that we expected classically. The nature 
of the two quantum phases other than the ferrimagnetic phase is not clear as of 
now. 

• For any 5 if all the spin-i interactions are made zero(J2 = 0) then the ground 

state energy must essentially be that of a spin-1 Hciscnbcrg antiferromagnet. In 
the thermodynamic limit this energy per site has been calculated [13] to very good 
accuracy. Our result {Eq/N = — 1.41712Ji) for the 8 site cluster agrees with the 
above value {Eq/N = — 1.40148Ji) upto finite size effects. 

2.3 Spin correlations in the ground state 

All the spin correlation calculations have been done with a system size of 10 triangles. In 
order to check the accuracy of the data the following checks were employed. 

• The three correlations < S^Sj >, < Si^Sj~ >,< Si~Sj~^ > were calculated sepa- 
rately. Then for each case it was checked that the latter two were equal and for 5's 
for which the cigenstatc was a singlet it was checked that the < ^j.^j > (which 
can be computed from the above three was thrice the < S^Sj > correlation as would 
be expected for a singlet. 

• For the singlet states it was verified that J2j < j > where j runs over all the 
spins for a given i was zero. 

In the graphs that follow, the < S^Sj > correlations are reported. The numbering 
scheme for the spins is same as shown in the figure in chapter 1. 
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2.3.1 Spin-1 - Spin-1 correlations in the ground state 



Spin-1 - Spin-1 correlations (S|^-S.^) in the ground state 
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2.3.2 Spin-1 - Spin-^ correlations in the ground state 



Spin-1 - Spin-1/2 correlations <Si^-S,^> in the ground state (l-spin1 ,j-spln1/2) 
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2.3.3 Spin-^ - Spin-^ correlations in the ground state 



Spin-1/2 - Spin-1/2 correlations (S.^-S.^) in the ground state 

0.06 1 , , 1 1 , , , , 




-0.1 2 1 ' ' 1 1 ' ' ' ' 1 1 

0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

5 

Salient features of the correlation function plots are the following: 

• Though for most of the region the correlation functions behave regularly, in the 
neighbourhood of 5 = 0.5 and d = 1.0 the correlations vary very rapidly with 6. 
From a preliminary examination of the data it seems that the rate of change of 
correlations with S is divergent at 5 = 1.0 . 

• The spin-1 - spin-1 correlations show the expected behaviour as we go towards the 
large 5 limit. The spin correlations are alternating in sign as we move from one spin 
to the next (going farther from the first spin) and also decaying with distance as we 
would expect for a pure spin-1 system. 

• The spin-1 - spin| correlations all tend to zero in the large S limit as we see from the 
plot. This is to be expected as in the large S hmit the system behaves essentially as a 
spin-1 system which forms a singlet and all the orderings of spin halves are essentially 
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degenerate as their interaction is much weaker in comparison to the spin-1 - spin-1 
interaction. Thus the spin-l's and spin-^'s are not strongly correlated and thus the 
correlations tend to zero. 

• The most interesting behaviour in the large S limit is shown by the spin-| - spin- 
I correlations. The most striking feature in this plot is that the next-nearest- 
neighbour correlation is larger than all the other correlations, even the nearest 
neighbour one. This and some of the other features can be explained if we cal- 
culate the effective Hamiltonian governing the spin-|'s in the large S limit. This 
calculation and the results that come forth from that analysis have been presented 
in chapter 3. 
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Chapter 3 



Perturbation theory and the 
effective Hamiltonian 

The model Hamiltonian we have been concerned with is the following: 

H — Jl^2 ln-~^ l{n+l) + J2 X!("^2n-'^ln + '^2n-"^l(n+l)) (3-1) 

n n 

where 5*1 's are the spin-1 sites and 5*2 are the spin-i sites. 

In the regime where the interaction between spin-1 's is much stronger than the in- 
teraction between the spin-|'s, we can consider the second term in the above Hamiltonian 
to be a perturbation to the spin-1 system. By doing perturbation theory calculations 
we can then find the effective Hamiltonian (for a few low lying states) governing the 
spin-i's once the spin-l's are essentially decoupled from them as a pure spin-1 system. 

We have found that to second order in perturbation theory this effective Hamiltonian 
has a particularly simple form which has only two spin interactions involving terms of the 
form In this chapter we set up the formalism to find that Hamiltonian and give a 

proof of the fact that to second order it has the form mentioned above. Then we describe 
a numerical technique to calculate the various interaction strengths in the problem. 
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3.1 The calculation of effective Hamiltonian 

We calculate the effective Hamiltonian in the following manner: 

• The first term of the Hamiltonian in equation 3.1 is treated as the unperturbed 
Hamiltonian Hq. And the second term is the perturbation V. Let > be the 
eigenstates of Hq, the spin-1 system. We assume that iV'i > are simultaneous eigen- 
states of Ho, the total angular momentum ((Zii^ii)^) the total Sz{J2iSli) 
operators. Such states can always be found as they form a mutually commuting set 
of operators. 

• The ground state of Hq is known to be a singlet and furthermore it is non-degenerate. 
We calculate the "corrections" to this ground state energy using non degenerate 
perturbation theory. 

• The perturbation term V contains both spin-1 and spin-| operators. But in calcu- 
lating the corrections to the ground state energy, the required matrix elements will 
be evaluated using the spin-1 system eigenstates. Thus we will be left at every order 
in perturbation theory with spin-| operators and their products. Just as we would 
have got the perturbative corrections to the unperturbed energy eigenvalues of Hq if 
the above mentioned matrix elements had been numbers, here we get perturbative 
corrections to the spin-1 Hamiltonian. These corrections order by order will consti- 
tute the effective Hamiltonian of the spin-| system. We note here that this effective 
Hamiltonian can be used to find only the states of the full Hamiltonian which lie 
close to the singlet ground state. That is because it is calculated by evaluating the 
effect of the perturbation only on the singlet ground state. 

Having set up the formalism we now proceed to calculate the effective Hamiltonian. 
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3.1.1 The effective Hamiltonian: First order 

The effective Hamiltonian to tlie first order (in J2 ) is given by, 

Ai/i = J2 < 0|(^( 5^2n ■ 5^1n + ~^2n ' ~S In+lW > (3-2) 
n 

wfiere we denote by |0 > the singlet ground state of the spin-1 system. 

Clearly, this is zero. That is because the state |0 > is a spherically symmetric state 
and the spin-1 operators occur in the above expression linearly. This can also be argued 
from Wigner-Eckart theorem. We know that all the spin-1 operators can be expressed as 
linear combinations of spherical tensors of rank 1. But the state with respect to which 
the expectation value is being taken is a singlet. Since we cannot add Ji — and J2 = 1 
to give Jtotai — 0, the above expression must be zero. 



So to first order in J2, tiie effective Hamiltonian vanishes. 



3.1.2 The effective Hamiltonian: Second order 

The effective Hamiltonian to second order in J2 is given by, 

where the A; 7^ implies that the sum is taken over all eigenstates except the singlet 
ground state. Ei is the energy of the state ipi and Eq is the energy of the singlet ground 
state. 
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Let us consider the matrix element, 



~^2n ■ ~^ln + ~^2n ' "^l(n+l)) | V'fe > 



The state on the left is a singlet (J = 0) and components of both and ~^i{n+i) can 
be expressed as spherical tensors of rank 1. Thus Wigner-Eckart theorem guarantees that 
the only those states \ipk > will contribute which have a value of total angular momentum 
which when added to J = 1 can give us Jfot — 0. But that means that the only allowed 
value is 1. Thus we come to the conclusion that in Eq. 3.3 we only have to sum over 
such > which eire spin 1 states. Equipped with this simphfication we now look at 
one particular term in eq 3.3 corresponding to a particular spin 1 state. It will look like, 



Eo—Ek 



where k labels the particular spin 1 state and the sum is over i which labels the particular 
Sz component (i=-l ,0, 1). 

One generic term in the above sum will look like 

^2 Eo-Ek 

We note here the the denominator will be the same for all the i's as the states have the 
same total angular momentum. Furthermore, we as of now don't make any assumptions 
as to the relative values of n and n'. Whatever we derive below will be true whether they 
are equal or not. In terms of components the numerator of the above expression will look 
like, 

where a , (3 = x, j, z. And < 5'"„ >q,. =< 0|S'"„|Vj > < >io is the complex 
conjugate of the same (for now we drop the superscript k as we will talk about a particular 
spin-1 state). 
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We make the following two claims: 



Ei < 5f„ 5f„, (A) 

• Ei < >oi< -^L' >io = if « # /3 (B) 



Before proceeding to prove the above we define the following: 

V2U+^ = -Si, = ^/2C/_l = S^^, x/2y_i = S^^„ U, ^ 5f„, ^ = S^^, 

U±ifl and V±i^o ^"^^ by definition components of spherical tensors of rank 1 and S'"„ and 
S'l^i can be expressed as linear combinations of components of U and V defined above. 
We now prove the above two assertions: 



Proof of (A) 

We consider the x-x term first. 

i 

+ < o| -^+1+^-1 >< ^-i r^+^^-^ O > (3.4) 
where we have summed over the three values of 5"*"* for the spin-1 state. 
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We can now use the m-selection rule to eliminate those terms above which are zero. 
After doing that we see that the above expression reduces to: 

i 

+ <0|^|^_i><^_i|^|0> (3.5) 
In an exactly analogous manner the y-y term reduces to: 

j:<si >o,< si, = - < oi^iV'i X Vil^lo > 

i 

-<0|^|V^_i><V^_i|^|0> (3.6) 
which is the same as the x-x term. Finally the z-z term is given by, 

E < >oi< SL' >io =< >< HVo\0 > (3.7) 

i 

Once we have proven that the x-x and y-y terms are equal the above has to be equal to 
the other two by rotational invariance. This can also be seen explicitly by application of 
the Wigner-Eckart theorem. The arguments are the following: 

• The matrix clement of the form < t/'jIOIO > wherever it occurs must have the same 
value everywhere as all the relevant Clebsch- Gordon coefficients are 1 (we are adding 
Ji = and J2 = 1) and the other term that we need to evaluate is the same for all 
such elements as it does not depend on the values. 

• The matrix element of the form < O|0|^i > even though has the same magnitude 
in all the three equations has the opposite sign in the z-z term, this is because 

< 1, ±1; 1, =f1|1, 1; 0, >= - < 1, 0; 1, 0|1, 1; 0, > (we have used the notation 

< Ji? Jiz! J2, j2z|ji, i2, i*"*, j*"* >)• This negates the sign difference in the right hand 
sides of Eqs 3.5, 3.6 and 3.7. Moreover the factor of | in the x-x and y-y terms is 
also accounted for by the fact the the x-x and the y-y terms contain the sum of two 
terms each of which arc equal in magnitude to the z-z term. Thus having proved 
that the x-x, y-y and the z-z terms are equal we have proved (A). 
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Proof of (B) 



We first consider the x-y term. We have after ehminating terms using the m selection 
rule, 

i:<SL >oi< Sl> >,o = + < Ol^l^i X ^il^lO > 

i 

-<0|^|^_i><^_i|^|0> (3.8) 

we see that both terms again have the same magnitude (note the point about the relevant 
C-G coefficients in the proof of (A)) but opposite sign so this term vanishes. 

Thus y-z and the x-z term must also vanish by rotational invariance (these can be 
also be shown explicitly using arguments similar to those used in the proof of (A)) . 
Thus having proved the two assertions we now come to the conclusion that the effective 
Hamiltonian governing the almost decoupled spin-|'s (close to the ground state) to the 
second order in perturbation theory is given by two-spin interactions of the form 
The final form will thus look like, 

J2 

Heff = a + -^[Ci(^i.^2 + ^2'^3-"-) 
+02(^1.^3 + ^2'^4"0 

+03(^1.^4 + ^2.^5...)] (3-9) 

where a = Na^J-i + N^bo and Ci, C2, C3 etc (upto a factor of J|/Ji) are the coupling 
strengths between the nearest neighbours, next-nearest-neighbours and so on. The first 
term in a corresponds to the energy of the spin- 1 system and the second term comes from 
spin- 2 terms of the form 

We note here that the q's involve a sum over matrix elements connecting all the 
spin-1 excited states with the singlet ground state (eq 3.3). Analytically calculating this 
sum is difficult as we do not have the complete information of all such states in order to 
calculate the required matrix elements. Thus we calculated the coefficients numerically. 
The method used is described in the following section. 
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3.1.3 Calculation of interaction strengths between the spin-|'s 

The computation of the interactions strengths was done using 6{= J1/J2) — 10 and for a 
system size of 8 triangles. This ensures that the perturbative corrections are convergent 
as the matrix elements calculated in eq 3.3 arc of the order J2 and for convergence of the 
perturbation theory this must be greater than the difference in the energy between the 
ground state and the first excited state of the spin-1 system. That is known to be of the 
order Ji. 

We calculated the interaction strengths (J|/Ji)ci, (J|/Ji)c2 etc as follows: 



• We begin by reducing all the bond strengths to the spin-|s to zero. This will give 
us the constant Oq. 

• Now we connect the bonds with strength J2 to only two of the spin-i s. We do this 
successively for nearest neighours next nearest neighbours and so on. 

• There being only two spin-i s in the system, the effective Hamiltonian governing 
them to second order in perturbation will be of the form A + 

• and being spin-| s we know that the ground state and the first excited states 
will have values A — ^B and A + ^B or vice versa. Which of these is the ground state 
depends on the sign of the coupling. If the coupling is ferromagnetic, the ground 
state will be a triplet and the latter will be the ground state and else the former 
will be the ground state. 

• Thus knowing the ground state and the first excited state and thus A and B we get 
bo and the q s. 



The various coefficients of eq 3.9 calculated using the above method (for a system size of 
8 triangles) turn out to be. 



ao = -1.41712, 


bo = 


-0.12665 


ci = 0.0183, 


C2 = 


0.1291 


C3 = -0.0108, 


C4 = 


: 0.0942 
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The important feature that we notice is that the next nearest neighbour cou- 
pling C2 is stronger than the nearest neighbour couphng Ci . This is the reason 
why the spin-| - spin-i correlations between the next- nearest- neighbours is larger than 
the correlations between the nearest neighbours. 

To this order in perturbation theory, the effective Hamiltonian seems to have the 
pattern AAFA (A antiferromagnetic, F Ferromagnetic). One curious thing that we 
notice in the spin-i - spin-| correlation plot is that though C2 and C4 are both positive 
(antiferromagnetic) the relevant correlations seem to be opposite in sign. The reason for 
this is that though both couplings are antiferromagnetic, C2 larger and thus it exercises 
greater control over the alignment of 6^5 than C4. That is why though C4 would dictate 
5^1 and 6^5 to be oppositely aligned, C2 being larger will win over in trying to align "^1 
and 5^5 in the same direction. This is what results in the correlations of opposite sign. 

For a system size of 10 triangles, the spins {1,3,5,7,9} and {2,4,6,8,10} because of 
a large C2 will form two distinct frustrated systems and thus the correlations reported 
for this size will be smaller than for system sizes which have even number of triangles 
(this has been checked for 8 triangles). The general features of the effective Hamiltonian 
nevertheless remain unchanged as we can see from the plots. 

Higher orders in perturbations theory: It is an interesting question to ask if the 
higher orders in perturbation theory will contribute in this case {5 — 10). Analytically the 
effective Hamiltonian will have a more complicated form than just the two-spin interaction 
and we cannot use the method just described used to calculate Cj numerically. We can 
try to eliminate the higher order effects by going to a higher 6, say 100. But calculating 
the eigenvalues corresponding to 5 = 100 to compute q to required accuracy is difficult. 
But if we look at the graph we see that even for 5 = 2 we already see the features that we 
have predicted for 6 = 10 quite clearly. This gives us an indication of the fact that the 
higher orders in all probability wont change the general features at 5 = 10. 
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Chapter 4 

Spin wave analysis 



In this chapter we look at the spin wave spectra obtained assuming the classical ground 
states described in Chapter 1 to be good approximations to the quantum ground states. 

4.1 Spin wave spectra in the Ferrimagnetic state 

The ferrimagnetic state which is stable for J2S2 > 2JiSi, is characterised by a collinear 
arrangement of spins with the spin-|'s and the spin-l's pointing along opposite directions. 
Thus the spin-l's and spin-|'s belong to two differnt sublattices having net magnetic 
moments in opposite directions. We proceed as follows to do the spin wave calculation 
for this phase, 

We write the Hamiltonian in the form, 

H^Hi + H2 (4.1) 

where Hi = JiY,n ~^ln ■ ~^lin+l), and H2 = J2 J2nC^2n ■ ~^ln + ~^2n ' ~^l{n+l))- 
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We now introduce the bosonic variables (we assume that the spin-l's are pointing 
along the +z direction), 

^i" = ^,al = \/25^&n (4.2) 

^In — Si — al^Qn, — —S2 + &j^&n 

which are essentially the Holstein-Primakov transformations for the two sublattices trun- 
cated to the lowest order (the equations for though are not approximations). And 
being bosonic variables [a„, al,] = 5„,„, and bl,] = Sn,nf- 

We also introduce the bosonic variables in the reciprocal space, 

«n = ^ E ^n = ^i: 6.e-^^("+^) (4.3) 

viV „. ViV „. 



We assume unit lattice spacing. Thus the wave vectors {k} go from — tt to tt in units of 
27r/N where N is the number of sites on each sublattice. 

Having defined these variables we now write Hi and H2 in terms of these bosonic 
operators assuming the classical configuration of the spins in the ferrimagnetic state to 
be a good approximation to the quantum ground state. If we look at Hi we see that the 
configuration of spin- Is mimics the ferromagnetic ground state of a Heisenberg ferromag- 
net in one dimension. So we can directly put down the form of Hi in terms of the bosonic 
operators taking into account the fact that in the present case Ji > 0. It will be, 

Hi = NJiS!-UiSiJ2sm''{^)alak (4.4) 

k ^ 

Now using Eqs 4.2 and 4.3 we write H^ too in terms of the fourier variables. This will 
have the form, 

H2 = -2NJ2S1S2 + 2J2S1 Ek blbk + 2J2S2 Efe 4«fe 

+2 J2^/;?^[Efc(&l4 + b^ak) cos f ] (4.5) 

Finally adding up Hi and H2 we get the full Hamiltonian in the fourier variables to be, 

H = NJiSl - 2NJ2S1S2 

+2J2S1 Efc blbk + Efc(2 J2^2 - 4 Ji5i sin2 |)4afc 

+ J:k{^J2^/S;S~2Cos miai + bkak) (4.6) 
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We note that the first term is the energy of the classical ground state. For convenience 
of expression we now define, 

= 272^2 -4Ji5isin2(^), Sfc = 272^1, Cfe = 2J2V^cos^ (4.7) 

{Bk is 2J25'i for all k). Thus we have, 

H = NJ.Sf - 2NJ2S1S2 + E Akaiak + ^ Bkblbk + E Ck{bWk + hat) (4.8) 

k k k 

We now use the Bogoliubov transformation to bring the above Hamiltonian to a more 
conveniently diagonalisable form. They are, 

flfc = Uk cosh 6k — v\ sinh 6^ 

bk — Vk cosh 6k — ul sinh 6k (4.9) 

Uk and Vk satisfy the same bosonic commutation relations as a and b. We for now leave 
6k as undetermined. Writing the Hamiltonian (4.8) in terms of the new variables Uk and 

Vk: 

H = -JrNSl - - Ek Ck sinh 26k + E(^fe + Bk) sinh^ 6k 

+ Efc(^fc cosh^ 6k + Bk sinh^ 6k - Ck sinh 26k)u\uk 
+ Efc(^fc sinh^ 6k + Bk cosh^ 6k - Ck sinh 26k)vlvk 

+ Efe[-(f + f ) sinh 26k + Ck cosh 2^^] {ukVk + 44) (4-10) 

The last term suggests that we choose 6k according to the following definition: 

IC 

tanh2gfc^ J (4.11) 

This reduces the Hamiltonian apart from constant numbers (which we can ignore by 
setting the zero of energy appropriately) to a form in which all the operators are of the 
form 6to. 

Thus the Hamiltonian finally looks like, 

H = E(Afc cosh^ 6k + Bk sinh^ 6k — Ck sinh 26k)u\uk 

k 

+ J2{^k sinh^ 6k + Bk cosh^ 6k - Ck sinh 26k)vlvk (4.12) 

k 

30 



where we have dropped the constants. So we have two modes given by, 
LO- = Ak cosh^ 6k + Bk sinh^ 9k — Ck sinh 29k, 

uJ-^- = Ak sinh^ 9k + Bk cosh^ 9k — Ck sinh 29k (4-13) 

The mode a;_ is gapless as k ^ and is also dispersionless when J2S2 — 2JiSi. The 
mode a;+ is always gapped as shown in the figure below: 



Spin wave modes in the ferrimagnetic state (J^=0.2, 0^=1) 




-3-2-10123 

k 



4.2 Spin wave spectra in the Spiral state 

We now investigate the spin wave spectra in the spiral state. This phase is stable for 
725*2 < 2/15*1. We assume a coplanar configuration of spins in the system as the classical 
ground state as shown in the figure of the same in Chapter 1. The basic unit of this 
ground state ( which keeps repeating) schematically looks as follows: 
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where cos^ = 72-52/2 Ji^*! and S\ and 5"^ denote spin-l's and spin-|'s respectively and 
the subscripts are the site indices. In order to do vahd spin-wave calculations with this 
classical configuration, for the spin-Is on the baseline we must choose(for each spin-1 
separately) the -|-z direction to be along (or opposite to) the classically expected direction 
of the spin vector. In doing so, we obtain the following transformation equations for the 
components of spin vectors (and thus for the corresponding operators) corresponding to 
the spin- Is, 



*-'2m 


= 52t'cos^-52'^sin^. 


qlx 

*-'2m-l — 


^2m-l COS 6 + S\^_ 


I sin 9 




q\z 


= 52t'sm^ + 52'^cos^, 

Q\y nVyl 
*-'2m — >->1mi 


qlz 

^2m-l ~ 


Qlx/ „• n 1 Qlzr 
'-'2m-l "T "^2™- 

qly qly/ 
^2m-l — ^2m-l 


_^ cos 9 


(4.14) 



Here we have chosen the direction in which spin-^'s are aligned to be the +z direction 
and the -|-y direction is into the page. The primed coordinates are obtained by rotation 
about the y-axis by an angle 9. Because of the different classical orientations of the spin- 
l's on the even and odd numbered sites we see that the transformation equations for them 
are different {9 has changed sign). We also note that the classical spin-1 vectors are along 
the —zf directions at each site. We can now proceed with the spin wave calculation as 
the approximation used (that the deviation of Sz from is small) is valid if we use the 
primed coordinates for the spin-l's. 



32 



As opposed to the ferrimagnetic case here we have to define three categories of 
bosonic variables. They are 





q1+ 
'-'2m- 1 


= ^j2Sid2m-l-, 




= ^2^24 






r.1- 
^2m-l 


= ^/2Sid2jn_i, 


r,2- 








QZl 

*-'2m-l — 


-Si + 0?2m-1^2m- 


c22 

-1) 


^82- clcn 


(4.15) 



Where the 6's d's are the bosonic variables for the spin-l's at the even and the odd 
numbered sites respectively (in the above expressions, etc have been defined in terms 
of the primed components). The index n runs over all the spin-| sites. Prom here the 
calculation proceeds in the following stages: 

1. We use Eq. 4.14 and Eq. 4.15 to write the Hamiltonian (4.1) in terms of the three 
bosonic variables in the unprimed coordinates. 

2. The first order terms in 62m's, c?2m-i's and c„'s vanish. We thus keep the terms to 
the second order in the above variables as the approximation to the Hamiltonian. 

3. At this point the Hamiltonian doesn't have a simple form which can be diagonalised 
using the Bogoliubov transformation used in the ferrimagnetic case. We use another 
method to find the spectrum here. We define a new set of canonically conjugate 
variables using the following equations, 

V2b2m = qb2m + ^PWm, V2&2m = Qb2m - m2m 

\Pid<im-\ = Qd{2m-1) + Wd{2m-1), V^4m-1 = Qd{2m-1) " Wd{2m-1) 

V2Cn = Qcn + ipcn, V^cl = ^cn " iPcn (4.16) 

4. Now we write the Hamiltonian obtained in step 2 in terms of these operators. In 
terms of these operators we find a couple of properties of the Hamiltonian. The q's 
and the p's don't couple in any of the terms. Secondly qb2m'^ and qd(2m-i)^ and the 
corresponding momenta occur completely symmetrically in the Hamiltonian. There 
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is no way to choose one over the other. So instead of having two different variables 
for the odd and even numbered spin-1 sites we can express the Hamiltonian using 
just one set of variables for all spin-1 's. We call that set Qin and Pi„ where n runs 
over all spin-l's. For convenience of expression we now call qcn,Pcn Q2n ^-nd 
respectively and n here runs over all spin-^'s. 

5. After all the above simplifications the Hamiltonian finally looks like (omitting the 

constant term), 

En A(Ql^ + Pi) + En 5(QL + PL) 

+ E„C'C0S26I Ql„Ql(n+l) + E„^C0S6' Q2n{Qln + Ql{n+1)) 

+ En C PlnPl{n+l) ' En D P2n{Pln + A(n+1)) (4-17) 

where, A — J2S2COS9 — JiSiCOs29 ~ JiSi (using the definition of cos^), B — 
J2S1 cos e,C^ JiSi and D = J2 ^5^. 



Clearly this is the Hamiltonian for coupled linear harmonic oscillators with nearest neigh- 
bour and next-to-ncarcst neighbour interactions. This tells us that for the purpose of 
finding the cigenfrequcncics we can use the classical equations of motion. This is because 
the eigenfrequencies of any such system on quantisation turn out to be the same as the 
classical ones. Thus using the Hamilton's equations of motion and the trial solutions, 

Qin = ^iq exp i{kn — out), Pin = eip exp i{kn — ujt) 

Q2n = ^2q exp i{kn - ujt) , P2n = ^2p exp i{kn- ut) (4-18) 

(where k e {— 7r,7r}) we get the following matrix equations, 

P„ = -A'Qn 

Qn = B'P„ (4.19) 



where. 



Pin 
P2n 



and Qn — 



Qin 
Q2n 



A' and B' are calculated to be. 
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A' 



2A + 2Ccos2^cos/c 2De-^ cos^cos f 



2L)e2 cos ^ cos I 



2B 



B' 



2^ + 2CcosA; 2De-^cosf 



2£)e2 cos I 



2B 



Clearly the eigenvalues of the matrix B'A' will give us the squares of the eigenfre- 
quencies. The following can be easily shown using the definitions of A,B etc. 

DetB'^0, DetA'^0 (4.20) 

Both of these together clearly imply that one of the eigenvalues of B'A' will be 0. 
Since the trace of B'A' will be the sum of its eigenvalues, the other mode can be calculated 
by calculating the trace of B'A' as one of the eigenvalues has already been determined to 
be 0. After taking the trace of B'A', we find the other mode to be given by 

(4.21) 



a;' = U^[{SiC0s9 - S2C0s''A) + &Un''0sm'' k]. 



The above mode has a minimum at A; = (a;(/c = 0) = 2 J2<S'2 1 1 — ^ | ) which vanishes 



at ^ = 2. At this point the nonvanishing mode looks as shown in the figure: 
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Thus the spiral state spin wave spectra consists of two modes one of which is dis- 
persionless and gapless and the other gapped with the gap vanishing at 5 = 0.5. Notably 
S — 0.5 is indeed one of the points where the system has been observed to be gapless 
numerically. Spin wave analysis doesn't give an indication of the point S — 1.0 being 
gapless. 

4.3 Conclusions 

To conclude, we have studied numerically and analytically, a two-spin (1 and 1/2) variant 
of the antiferromagnetic sawtooth lattice. Interesting features brought to light by this 
study are the following: 

• The system seems to be gapless at two points in the phase diagram. Thus there seem 

to be more than the two expected phases from the classical analysis. Moreover none 
of these points correspond to the classically expected point for the transition from 
the ferrimagnetic to the spiral phase. Spin wave analysis does give some indication 
of the point 6 = 0.25 being gapless but the points which are actually found to be 
gapless are S = 0.5 and 6 = 1.0. This may be due to the small values of spins 
because of which spin wave analysis may not be very accurate. The properties of 
the system at 5 = 0.5 and 5 = 1.0 are unclear as of now. 

• In the large S limit, though the other correlations behave as expected, the corre- 
lations between the spin-|'s have the curious feature that the strongest coupled 
spins are the next-nearest-neighbours. This we have been able to explain using 
perturbation theory analysis of the model. 

The study of the nature of the quantum phases at the points S — 0.5 and S — 1.0 
is an interesting direction in which further work on this model can proceed. Numerically 

one can try to go the larger system sizes using DMRG to eliminate the finite size effects 
and thus better approximate the thermodynamic limit. The aspect of the problem we 
haven't touched on at all is the effect of a magnetic field and the thermodynamics of the 
system. This is a another direction of study which can lead to a better understanding of 
this model. 
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